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We extend a perturbative, nonequilibrium renormalization group approach to multi-orbital sys- 
tems and apply it for studying transport through two parallel quantum dots coupled electrostatically. 
In general, the conductance shows pronounced Kondoesque peaks at three voltages. One of these 
peaks disappears if, as in some experiments, one of the dots is decoupled from one of the two leads. 
For an asymmetric coupling to the leads, also negative differential conductances are possible. This 
is a genuine nonequilibrium effect, accompanying the Kondoesque peaks. Moreover, a criterion to 
distinguish spin and orbital Kondo effect in such a system is discussed. 



PACS numbers: 73.63.Kv,71.27.+a,72.10.Fk 



I. INTRODUCTION 

Kondo physics in quantum dots has been under in- 
tense investigation in recent years. As was predicted 
theoretically^^ the Kondo effect can lift the Coulomb 
blockade at low temperatures, leading to a conductance 
peak at zero bias^ In such experiments, often more than 
one orbital is involved. This is unavoidable, if the tem- 
perature or Kondo temperature is larger or comparable 
to the mean orbital splitting. Some nanostructures are 
also intendedly constructed in such a way, for example, 
two parallel quantum dots coupled electrostatically-' 5 ' 6 
The latter allows for a controlled study of the interplay 
of spin and orbital (upper/lower dot) degrees of freedom. 

On the theoretical side, the TV-fold degenerate situa- 
tion where all N (spin and orbital) levels of the Ander- 
son model have the same one-particle energy e, Coulomb 
interaction U, and hopping amplitude t to the leads is 
best understood. At strong coupling, the Schrieffer- Wolff 
transformation 7 then leads to the Coqblin-Schrieffer 
model 8 in terms of the exchange coupling J. This model 
has SU(N) symmetry for the rotation of orbitals and spins 
and, with increasing degeneracy N, leads to a strongly 
enhanced Kondo temperature j^ u W w j^ e -i/(Np j) 
(Pq: density of states in the leads; D: bandwidth of the 
leads' we set ks = h = e = 1 unless these constants 
appear explicitly in the equations). 

Of course, we cannot expect a real quantum dot 
to have SU(N) symmetry, unless special arrangements 
are madeiii Hence, quite an effort was devoted in the 
Iiteraturei2ii&i2ii£ii& to study the situation where ei- 
ther a magnetic field or a difference in orbital energy 
8 splits N = 4 levels into two doubly-degenerate sub- 



sets, hence breaking the SU(4) symmetry. If 8 <C T^ U( - 4 - ) 
the low-energy physics still resembles that of the SU(4)- 
symmetric model. In contrast for 8 ;g> one of the 

levels drops out of the scaling procedure; the low-energy 
physics is that of the usual SU(2)-symmetric Kondo 
model. In between, there is a crossover region^ 5 - 

It is also unavoidable that the coupling constants J are 
orbital-dependent in experiment. Numerical renormal- 
ization group (NRG) calculationsiiii 5 " suggest that for 



small enough 8 both orbitals are screened at the same 
energy, in agreement with scaling analyses showing a ro- 
bust (marginal) SU(4) fix point. 15 ' 16 ! 17 

Let us also distinguish here between orbital conserving 
and non-conserving dot-lead couplings. In our paper, we 
consider the former, experimentally realized by separate 
leads for each dot^ see Figure The non-conserving 
case is obtained if the quantum dots couple to a single 
lead like in Ref. 0- Then the hopping processes do not 
conserve the orbital quantum number, which can give rise 
to SU(2) Kondo physics even if more than one degenerate 
level is involved, see Ref. ^| and references therein. 

So far, we only discussed equilibrium physics. A fi- 
nite voltage V leads however to genuine nonequilibrium 
effects. For example, if V » 8 ^> Tk (the Kondo temper- 
ature) scattering between the two subsets of levels be- 
comes relevant again. A Kondoesque resonance develops 
at V~8, visible e.g. as a peak in the conductance. But, 
at the same time, a genuine effect of nonequilibrium is the 
presence of decoherence processes induced by the finite 
current. This decoherence cuts off the flow to the strong 
coupling fix point even for temperatures T <IC Tr- 
Then, physical quantities such as the differential conduc- 
tance depend strongly on the configuration of the system 
like the asymmetry of the hopping t between quantum 
dot and leads. 

Nonequilibrium Kondo physics is not as thoroughly 




FIG. 1: Sketch of the parallel double quantum dot config- 
uration. Each dot is coupled to its own leads; the coupling 
between the the dots is only via the Coulomb interaction (not 
indicated in the Figure). 
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investigated theoretically because standard approaches 
for the Anderson impurity model like NRG are not 
applicable. Rosch and coworkers^ recently developed 
a perturbative, noncquilibrium renormalization group 
(RG) method to address nonequilibrium Kondoesque 
physics^ This approach is similar to Anderson's^ poor 
man's scaling approach but includes nonequilibrium ef- 
fects caused by the finite current through the system. 
In nonequilibrium, the low energy physics cannot be ab- 
sorbed properly by a pure renormalization of the bare 
couplings. It is necessary to include the frequency de- 
pendence of the vertex functions. 

The outline of the paper is as following: In Section 
ITT1 we extend the perturbative RG of Rosch et a/iSi to 
many orbitals, supplemented by the calculation of the 
decoherence rates in the Appendix. In Section II I II we 
apply this method to the parallel double quantum dot 
system, with each dot being coupled to separate leads. 
Thereby, the first part, Section IIII Al discusses briefly 
the simplification of the RG Eqs. for this special applica- 
tion. The following Sections present our results: Section 
IIII Bl shows the renormalized vertex function, primarily 
of theoretical interest, Section Till CI the increase of the 
decoherence rates with voltage, and Section IIII Dl the 
change of the occupation in non-equilibrium. The reader 
who is mainly interested in the differential conductance 
may also directly turn to Sections IIII El and IIII Fl which 
present our results for this physical observable in the case 
of symmetric and asymmetric coupling to the two leads, 
respectively. In the case of symmetric couplings (Section 
IIII Efl we generally obtain two orbital Kondoesque peaks 
at V « ±<5 in the differential conductance of both dots 
and an additional spin Kondo peak at V = for the dot 
which is lower in energy. Depending on the parameters, 
the peak at V = can be strongly suppressed and hardly 
discernible. The main findings of Section UlI Fl for asym- 
metric couplings are: (i) in very asymmetric situations 
the Kondoesque peaks at V » —6 can disappear and (ii) 
negative differential conductances are possible for V > 6. 
In addition to these non-equilibrium results, the linear 
(equilibrium) conductance is analyzed in Section IIII Gl 
A summary of the results can be found in Section IIVI 



II. NONEQUILIBRIUM PERTURBATIVE RG 

Starting point for modeling a quantum dot with N- 
levels should be the Anderson impurity model: 

H AlM = H + ^2(t Xm c+ kX d m +H.c.) 

Xmk 

^ ^ ^ m^m^m ^ ^ ^ rnrn' T^rnT^rn' • (1) 
m m^m f 

Here and d m (c^ fcA and c m kx) are creation and an- 
nihilation operators for electrons in the dot (lead A); 
n m = d m d m ; e m describes the one-particle energy of level 
to with to = {cr, i} subsuming the spin and orbital index; 
Umm' is the Coulomb repulsion within the dot. The levels 
of the dot hybridize via t Xm with non-interacting leads A 
described by (later, we simplify e\ m k = £&) 

Ho = e \mk C \mk C \mk- (^) 
Xmk 

In this paper, we assume each quantum dot level to to 
couple to its own lead channel which is for example the 
case if the different levels correspond to different quan- 
tum dots like in Ref. (see also Figure ^ of the present 
paper as an illustration). But in other situations it is pos- 
sible to move an electron phase coherently from one level 
to via the leads to another level to'. This goes beyond 
Hamiltonian QJ. 

While the Anderson impurity model allows us to make 
contact with experiment (estimating parameters), we ad- 
dress in the following the parameter regime e m -C fJ-\ <C 
U m m' + &m f° r all to, to', A (fiy. chemical potential in lead 
A) so that charge fluctuations are suppressed. Then, with 
^ m (ti m ) = 1 electron within the dot, we can map Hamil- 
tonian ||TJ onto the subspace with one electron in the dot 
by a Schrieffer- Wolff transformation^ i.e., onto the gen- 
eral effective Kondo Hamiltonian (neglecting a potential 
scattering term) 



tt _ tt , Y-mm , / tAiA 2 vm 1 m 2 + , fAiA 2 v-miTOi_+ r \ f3l 

n — a a + 2.^ m + \ J mim 2 ^ C X 2 m 2 k2 C X 1 m 1 k 1 + J m im2 ^ C X 2 m 2 k 2 C X 1 m 2 k 1 > • W 

m X\m\k\^X 2 m 2 k 2 



Here, the Hubbard operators^ I™ 1 " 12 represent the 
scattering of the local state from level m-i to level m,\. 
When the lead's band edge D > \e m \, the local level e m 
of the Anderson impurity model is renormalized to e rn in 
the Kondo model; 10 J^ 2 = t\ im t x ^ ( Um ^ 2+£d ~ £) 

and J^ 2 = (1 -2S mi m 2 )t* Xim2 t X2m2 Umi * 2+ed where e d 
is the average level energy, neglecting the level splitting 



which is much smaller than U, €d- Note, that for mi = TO2 
both J and J yield the same kind of contribution. The 
advantage of having this term twice (i.e., distributed to 
mlwi an d ^mimj ls * ne simplification of the following 
equations. The same distribution of the toi = TO2 contri- 
bution is usually also employed for the Coqblin-Schrieffer 
model which is the SU(N) symmetric version of Hamil- 
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FIG. 2: One- loop Feynman diagrams for the renormaliza- 
tion of the vertex functions (a) ff m 1 1 m 2 (w) and (b) g^£ 2 (ui). 
The full line denotes the electron and the dashed the local- 
ized state of the Kondo model. Note that the second (hole) 
diagram of (b) does not contribute for (a) since it cancels 
with the corresponding first (particle) diagram. For the ver- 
tex renormalization, the electron line is restricted to a small 
interval \dD\ at the band edges, as indicated by the dash. 



tonian © [J and J = J/N are then level-independent]. 8 
Let us now generalize the perturbative RG approach 
of Ref. |^] to the multi-orbital Hamiltonian J3J . The base 
quantity is the vertex function which -at finite voltages or 
currents- becomes frequency dependent. Like in the poor 
man's scaling approach, we integrate out high energy de- 
grees of freedom in the intervals [—D,—D + \dD\] and 
[D — |cLD|, D], resulting in a renormalization of the ver- 
tex. The Feynman diagrams for this renormalization are 
shown in Figure to lowest order (one- loop RG). To di- 
rectly apply standard perturbation theory, the localized 
states of the Anderson impurity model can be expressed 
by pseudo Fermions, see Ref. [2l|. Starting point of the 
RG scheme is the unrenormalized ^-independent vertex 



"mi m 



M = PQ-Cima) Where We 

neglect the energy dependence of pq. Here, the frequency 
of the localized state is approximated by e m so that the 
vertex depends only on one ui (the average of incoming 
and outgoing electron frequency). This is only possible 
in the weak coupling regime. 

For the nonequilibrium situation, the calculation of the 
diagrams of Figure [21 requires in principle the usage of 
KeldyshSi Green functions. Rosch et alS^ argued how- 
ever that, to leading order in 1/ Ih(V/Tk), it is sufficient 
to keep track of the real part of ffm 1 1 m 2 (w), as indicated 
by perturbation theory^ This gives rise to the same kind 
of vertex renormalization per energy interval In D as in 
the poor man's scaling approach, i.e, at T = 0: 



1 8 f D , signje) 
2dlnDJ_ D V 



Auj 



(4) 



Here, Auj is given by the chemical potential and the 
change of energy of the localized state in Figure [21 A 
finite temperature, smears out the sharp step of the in- 
tegrand in Eq. cutting off the emerging logarithm. 
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FIG. 3: Vertex corrections and self energy diagrams beyond 
the one- loop approximation of Figure [21 
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FIG. 4: A similar combination of vertex corrections and self 
energy diagrams as in Figure [21 emerges if we calculate gener- 
alized spin/orbit susceptibilities. 



Another cut-off is induced by the finite voltage, driving 
the system out of equilibrium. This nonequilibrium ef- 
fect originates from higher order diagrams such as those 
of Figure [3 These diagrams give rise to an imaginary 
part of the self energy and vertex corrections, even at 
low w's for the nonequilibrium situation. Generally, this 
will result in the following modification of the pole in 
Eq. Q: 



1 d 



2 d In D 



D 

de- 

-d e 



sign(e) 



Auj + i7sign(e) 



7 2 ) 



Of 



(5) 



where 9(x) = 1(0) for x > 0(< 0). 

For the proper calculation of the low energy cut-off 7, 
we would need to calculate diagrams as in Figure^ This 
goes beyond the scope of the present paper. As in Ref. 
l2ll we follow a heuristic approach instead. To this end, 
first note that similar combinations of vertex corrections 
and self energy inclusions occur when calculating gener- 
alized susceptibilities for level mi and 7712; see Figure 0] 
Secondly, an incorrect prefactor of 7 is not affecting the 
results to leading order in l/ln(V /Tr-). Hence, we can 



crudely estimate 7 by the decohcrcnce rates 7„ 



of the 



corresponding susceptibilities. For example, we can take 
the maximal r y mi m 2 of the levels involved (or tempera- 
ture if higher): 7 = 7 mim2m3 = max{T, 7 mim2 , 7 m2m3 } 
where mi, m2 and UI3 denote the incoming, intermedi- 
ate, and final local state. In our calculation for two levels 
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in Section ITTT1 with SU(2) spin symmetry, only two of the 
orbitals i (m = {cr, i}) can be different, and the vertices 
are symmetric w.r.t. the spin indices. We then actually 
use the intra-orbital rate 7 = max{T, 7^ = 7^} as a cut- 
off for diagrams where only one orbital is involved and 
the inter-orbital 7 = max{T, 712} whenever both orbitals 



are involved. The calculation of r ) mi m 2 via the suscepti- 
bilities is presented in the Appendix. 

Having identified the cut-off r y mi m 2 m 3 , we can now 
write down the RG equation for the vertex functions 
ffm, A m 2 ( w ) and 9m?m 2 (u), corresponding to Figure 



^mat^) _ \^ AiA / , e ™i e "2 \ AA 2 / , e TOi £ W2 mT^m^a 

Q]jlD ym 2 m 1 \ LO ^ 2 /•yTci\Tti 2 \ ' ^ ' w+£mi — «m 2 — MA ' 

^ma^) _Y^n A i A r.- : £?n2 - £ "M „ AA 2 /, . e '»i ~ £ "t Np,7 m2mmi 
m A 

, Vn AlA C, iUn AAa fr) £ff "~ Em M rf m3TOl '"' « AAa 6 1 £ "i ~ £ "2 \q7^-2-i 1 

+ Z^ ff "H™2 l W iL3m im2 l w 9 i L L-(£ mi - Cm2 )/2- AtA _ 3m2m2 i W 2 J L L-(e mi -e m2 )/3- M J 

A 

l n AA 2 /,.\WiA I . £ mi -£m 2 N, p( 7m 2 m 2 m 1 ~Ai A / , e mi — £m 2 s p .7m 2 m 1 m 1 ] 



The equation for g can be obtained from Figure^- The 
first line of the equation for g corresponds to the right 
diagram of Figure[2b; the second and third line stem from 
the particle and the hole diagram of Figure^ with m = 
m 2 and m = mx, respectively. Note that the different 
frequencies of Figure|3are related by energy conservation 
and that ui = (wi+a;2)/22I. The fc-integration (note that 
k is not conserved) over the band edge states results in 
a factor po which is already included in g. The electron 
states in the leads range to ±D relative the the respective 
chemical potential fj,\. 

The renormalized vertex functions of Eq. Q have to be 
calculated self-consistently together with 7 mm < [Eq. I|A.7|) 
of the Appendix] which determine the cut-offs '}m 1 m 2 m 3 - 



Conductance 

To study the most important physical quantity, 
the conductance, we employ the scattering T-matrix 
T^; 1 i A ^ 2 (fci, fc 2 ) for the process moving an electron from 
state {fcimi} in lead Ai to state {fc2?7i2} in lead A 2 . This 
process necessarily also involves the local state which 
changes in the opposite direction. The change of the local 
state can be expressed via, spin and orbital operators in 
the T-matrix (see e.g. Ref. Eg). By (|f^ 2 (fci, fc 2 )| 2 ), we 
then denote the thermal average of the T-matrix w.r.t. 
the local state configurations. The steady state condition 
requires in terms of this averaged T-matrix: 

E< l f »n™2 ( fc i . *a) D/AifoJ [1 " Mtk 2 )] = 

Ai A 2 ;feife 2 

Ed^^.^jiv^yii-Afe)], (7) 

AiA 2 ;feife 2 



where f\(c) = 1/ '(e^ e_AlA ) + 1) is the Fermi function for 
lead A. A difference to Ref. I2lll30l is that we do not 
take into account a 7-broadening of the spectral func- 
tion. This is effectively descri bed by a broadening of the 
Fermi functions in Ref. 1211301 resulting in corrections to 
subleading order in 1/ \u(V/Tk)- 

Since we consider the situation where either the 
voltage-induced decoherence rate or the temperature is 
larger than the Kondo temperature the renormalization 
process is cut off, instead of flowing to strong coupling. 
Therefore, we are at weak coupling and we can calculate 
the scattering matrix in lowest order: (fci, k%) — 

(A2f7i2fc2|^int|Aimifci), where H m t denotes the interac- 
tion between the lead electrons and the local states via J 
and J, i.e., the last two terms of Eq. (J3J). Averaged over 
the local states, this yields: 

<|tf^(*i.*a)| 2 > = (1 - S mim2 )(n m2 )9 X m \tM 2 /p 2 o 

+S mim2 ( n ™)3mimH 2 /P0 
m/mi 

+S mim2 (n mi ) (g^X (uj) + g^ 2 H) 2 / >g. (8) 

Here, we have replaced the bare coupling by the renor- 
malized one, i.e., J — > g(uj)/po with u> = (e& + e^t)/2; 
(n m ) is the average occupation of local level m. 

In a first step g and g are calculated by Eq. 10, which 
is solved self-consistently together with the cut-off 7 mim2 
[Eq. (|A.7(l ]. With these g's and g's, the occupations (n m ) 
can be determined via Eqs. Q and JSJ and the constraint 

From the T-matrix [Eq. (jSJ)], we can then calculate 
the current from lead A and orbital m as the difference 
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between scattering out of and into the lead: 



Ixr, 



-27T <>(e k ~e k ' + e m '-e m ) 



Am' (^Am) 
kk' 



{/A(e fe )[l-/A'(e fc ')](&(fc,fc')| 2 } 
-h{e k/ )[l-h(e k )](\f^ m (k',k)\ 2 )}.(9) 

The differential conductance follows as G\ m = 
dhm/dV. 

III. TWO PARALLEL QUANTUM DOTS 

Let us now apply the perturbative RG equations to 
the situation of two parallel quantum dots each coupled 
to two leads, see Figure ^ of the introduction. For the 
Anderson impurity model, this dot-lead coupling is de- 
scribed by hopping processes txi] additionally there are 
intra-dot and inter-dot Coulomb repulsions Uu>. We 
study here the corresponding Kondo model which is re- 
stricted to the subspace with one electron in the dot. 
Then, g^, describes the combination of two hopping 
processes, leaving the number of electrons in the dot un- 
changed. 

A. Simplification of the RG Eq. 

The two parallel quantum dot case corresponds to 
quantum numbers m — {i,<?} with a —\ (+), j (— ) de- 
noting the spin and i = 1(+), 2(— ) the two dots (orbitals) 
with level energy = —iS/2 (i = —i in the following). 

Due to the SU(2) spin symmetry in each dot, we can 
reduce the number of vertices we have to deal with. Fur- 
thermore, some of the vertices are connected via scaling 
invariants. Altogether this allows for the following re- 
duction of vertices: = g* x 'H/2, g XX q„,(u) = 

{»)/*, SgtsM 

-XX' 



g"" (oj) /4, and -from the scaling 



invariants- g xx ! itT , (w) = [g xx {u)-g^ x (w))/4. The mean 



ing of the redefined vertices becomes clear if we write the 
Hamiltonian in terms of the coupling constants Jj(i2) and 
J corresponding to <?j(i2) and g: 



H = E( £fc ~ Ma) 



Xicrk 



c \iak c \iak ~ ^{ni~n 2 ) 



(10) 



jA'A 

El 7T~ • c\ iak T aa tC\ii a t k t 



XX' iaa' kk' 



E 



AA' 00' ' kk' 



tX'X 

2 



c A2 CT fe c A'l ff 'fc' + H.C. 



j 

+(ni-n a ) -^-( c xi<Tk c y^k'-ci 2ak cx'2ak')- 



A A 



XX' crkk' 



Here, Si denotes the spin operators in orbital i, t the 
vector of Pauli matrices, and n, = Eer^t^- Note that 



with only one electron in the dots, there is no need to 
consider Hund's coupling. 

Using the SU(2) symmetries in both orbitals and the 
above scaling invariants, the scaling Eq. 10 become par- 
ticularly simple for a symmetric coupling to the leads 
(then, g%((jj), gi(co), and gi 2 (co) do not depend on the 
lead indices; we can hence drop the A, A' indices): 

dgjjuj) 
dlnD 

dgiju) 
dlnD 

dlnD 



-E^( w ) :VA^/2+.9? 2 (^+f)e3^ (5+Ay/2) ] 



A=±l 



= -512M 



A=±l 



53 <&(w + ity2)e™ 



uj+i(8+XV/2)- 



A=±l 



(11) 

Eq. has to be solved self-consistently together with 
the decoherence rates of the Appendix which simplify for 
the parallel double dot system to: 

7< =7r£ fdu {^(u)/x(w)[l-A'(w)] 

AA' J 

+ &( u )f x ( u+i l)[l-f xl (u,-ii)]} (12) 



712 



XX' i 



IE. l<^{[-/i{") + ~9K")}fx{u)[l -/a'H] 



+g l 2 {u)fx{u+i~)[l-fx>{u-i~)]} (13) 
These 7$, or T for T > 7$, cut off the O-functions in Eq. 

CD- 

The equations for the occupations [Eq. Q, taking 
m\ = m 2 in Eq. and the current [Eq. ©] reduce 
to 

(m) £ I degl 2 (e-i-)fx(e)[l-fx>(e-i5)] = (i - i) (14) 

A, A' 

and 

I ia = f/dw{[3(n i ) ft 2 (w)+9?(a;)][/ i (a;)-/ Jl (a;)] (15) 
+4(r h )g 2 12 (u+i5/2)f L (u J )[l-f R (u>+i5)} 
-A{ni)gl 2 (u-i8/2)f R (u)[l-f L (u-i8)] 
+i(n,)g 2 12 (u;+iS/2)f L (u J )[l-f L (u J +iS)} 
-A(rH)&(u-i5/2)f L (u)[l-f L (u-iS)]}, 

respectively. In Eq. (|15|l . the first line stems from the 
current from left channel i to right channel i. The 2nd 
and 3rd (4th and 5th) line correspond to the current to 
right (left) channel i. The latter are balanced by similar 
terms from the right lead so that the net current from 
one orbital to the other is zero. 



B. Results: Renormalized vertex function 

Let us now discuss the numerical solution for the renor- 
malized vertex functions of Eq. JSJ / l|ll|) which is pre- 
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150 



FIG. 5: Renormalized vertex functions to the left lead at 
T/T% = 0.5, V/T% = 20, 8/T% = 65 {g l(12) in the legend box 
stands for g^n 2 ) )• Part a) is for symmetric initial couplings 

gi X ' = 0.04, g&' = 0.06 for all A, A'. In part b) the right 
lead of dot 2 was removed, i.e., parameters as in a) except for 
32(12) = 92(12) = 0- In part c) also the coupling of dot 1 is 
g{ L = 0.04, 5l M = 0.06, gh L = 0.04, g? R = 0, 



asymmetric: 
gi 2 L = 0.08. 



sented in Figure [5] This vertex function, a frequency- 
dependent renormalization of the interactions J, J be- 
tween lead electrons and local states in the quantum 
dot, is the fundamental theoretical quantity; it is not di- 
rectly observable experimentally, but -together with the 
occupations- allows to calculate e.g. the current which 
will be discussed later. In the following, all energies are 
normalized to the (equilibrium) Kondo temperature for 
6 = 0, i.e., to T£ = T K {5 = 0). 28 Due to the finite voltage 
V applied, the two leads are thermalized by two different 
Fermi functions with [il = fj, + V/2 and fiR = fi — V/2, 
respectively. The quantum dot is out of equilibrium. 

Here and in the following figures, the initial coupling 
constants which are not explicitly mentioned follow from 
the ones of the Figure caption since they stem from the 
same t's of the Anderson impurity. For example, g\^ = 
9i2 V9i R /9i L - We always set g* x ' = initially. 

To understand the four peak structure of Figure [5] it 



is instructive to integrate Eq. (|llfl analytically, dividing 
the Eqs. by the squared coupling constant from the left 
hand side (e.g. multiplying the first line by l/gf(uj)) and 
afterwards keeping the c/'s on the right hand side con- 
stant (at the initial values). This is justified as long as 
the renormalization of the vertex functions stays small. 
But our mere aim here is to understand the positions and 
widths of the different peaks. From Eq. (|llfl we straight- 
forwardly get in this way: 



912 M 
9i{u) 



1 



Y^Qjx l"+y 2 l' + in \_i__i_±_L^ 



I (XV-iS) | 



^ m \u+i{S+XV/2)\ 12 ' 



(16) 
(17) 

(18) 



where \u>\i = \j uP- + max 2 {7i, T}, and prefactors of or- 
der one (like gnl 9i) before the logarithmic terms are ne- 
glected. 

By means of Eq. I|16(l , we can now identify the origin 
of the four Kondoesque peaks in the g^s seen in Figure 
|5] When the frequency matches the difference to the 
left or right Fermi level (u> = ±V/2) a single dot (spin) 
Kondoesque resonance develops; but the flow to strong 
coupling is cut off by the corresponding decoherence rate 
7i (or T). The other two peaks in gi stem from the or- 
bital Kondo effect. Here, the localized state changes to 
the other orbital and back to the original orbital. These 
processes are enhanced if the conduction electron scat- 
ters to the Fermi level of the left or right lead, requiring 
u> = —i5 ± V/2. Hence the resonances [and those of Eq. 
(|18fl ] arc at ui = —id ± V/2 and the cut-off is 712. Simi- 
larly, the peaks for (w) in Eq. (|T7|) are at ±^-^ (with 
6/2 instead of S since u> is the average of incoming and 
outgoing frequency which differ by 5 as 512 changes the 
orbital) . 

In part b) of Figure [S] the right lead of dot 2 was re- 
moved. As scattering to the removed lead is now prohib- 
ited, one out of four peaks disappears for gi and g\2- The 
missing peak is the one at the lowest value of u> in part 
a). We can hence identify the missing peak with spin 
(52) and orbital scattering processes (g± and (712) via the 
right lead of dot 2. Part c) of Figure[S]is for the situation 
where also dot 1 is coupled asymmetrically. Some peaks 
are clearly suppressed. 



C. Results: Decoherence rates 

Figure El shows the three decoherence rates 71, 72, and 
712 which need to be calculated self-consistently together 
with the RG Eqs. As in Eqs. lfTS )l -l(T5 )l it is instruc- 
tive to calculate 7 from the unrenormalized g's (not self- 
consistently/in lowest order approximation). At large or- 
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bital splitting |<5| 3> max{T, this yields for symmet- 
ric coupling to left /right lead (Figure El around zero): 



\V\ • 



71 max{T, — }, 



72 
712 



Wl2\?>\, 
2^12^1- 



(19) 



The decoherence rate for dot 2, i.e., 72, and the intra- 
orbital rate 712 are both proportional to |<5| as there are 
many decoherence processes, reflecting the instability of 
the high energy level of dot 2. Decoherence processes for 
dot 1 on the other hand become only available at finite 
bias voltages. For very large \8\ 3> the system shows 
the usual single-orbital (spin) Kondo effect. Note that, 
for small V, we have 71 < T^. One might envisage that 
this signals a flow to strong coupling and, hence, a break- 
down of the RG Eqs. However, for the large splitting 5 
the real Kondo temperature is very much reduced. It 
actually is only m 1/100 here, much smaller than the 
cut-off. 

Let us now discuss the large voltage range of Figure 
El For \V\ » I 8 I, all curves show a similar asymptotic 
behavior: 

j i ^M(g i ) 2 + (g 12 ) 2 ]max{T,^-}, 

3 \V\ ( 2 °) 

7l2 « 2 7 r^[ i ( ffl ) 2 + (kf + (.9i2) 2 ] max{T, U}. 



In between the two extremes \V\ <C \8\ and \V\ ^> 
\S\, the transition from a strong orbital decoherence rate 
of dot 2 due to the high level energy to strong current 
induced decoherences in both dots is clearly visible at 
V « ±5 in Figure El 

Neglecting the finite orbital splitting, we can also get 
the higher order asymptotic behavior for the decoherence 
rates analogous to Ref. ITil 



7T 

7i ~ 7 



V 



4 (In 



v 

^ XT 



111 ■ 



(21) 



D. Results: Non-equilibrium occupation 

The numerical results for the different electron occupa- 
tion of the two orbitals is shown in Figure [7| for different 
orbital splittings 8. To ensure that the RG approach is 
valid, i.e., that we stay in the weak coupling regime, we 
use here and in the following a rather high temperature 
T/Tft = 2. Such temperatures are actually relevant in 
many experiments where often the lowest accessible tem- 
peratures are of the order of the Kondo temperature. At 
low bias voltage V < 8, only the lower orbital is occupied 
for sufficient large 8, as is to be expected for 6 ^> T. At 
high bias voltage V > 8, the electrons tend to stay in 
both dots with the same probability; (m — 71-2) decreases 
with increasing V. 



While all 7i's are now proportional to \V\, we see in Fig- 
ure El that nonetheless 7* <C \V\. This is due to the 
small prefactors ~ g 2 in Eq. (|20() . The current which 
has similar prefactors causes the decoherence processes. 
For sufficiently large V, all 7i's become larger than tem- 
perature; the current induced decoherence exceeds the 
temperature effect. 
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Scattering rates as a function of bias voltage V at 
0.2 and = 40. The initial couplings are sym- 

0.04, g$ x ' = 0.01, gii' = 0.06. For large bias 
voltage, the scattering rates are larger than T (dashed line) 
but much smaller than V. 



E. Results: Conductance for symmetric coupling 

Let us now turn to the physical observable of inter- 
est, the differential conductance through the two dots 
Gi = dIi(V)/dV (i = 1, 2), presented in Figure El The 
conductance consists of the intra-dot spin scattering and 
the inter-dot orbital scattering. For 8 — 0, a pronounced 
peak is found at V = 0, characteristic of the orbitally 
enhanced Kondo effect. The flow to strong coupling is, 
however, cut off by temperature. If the orbital degener- 
acy is lifted, the peak at V = is suppressed. For dot 
2, it finally disappears at large 8 when this dot is de- 
populated. For dot 1, we still have a spin Kondo effect 
at large <5, but with a very much reduced Kondo tem- 

hSU(2) 



perature Tj^ 2 -' 



« T°k- 



Because also T K u < T, the 



peak at V = is very much smeared out. At finite 8, 
two new peaks occur at V = ±8, corresponding to or- 
bital scattering. For T/T^ = 2, we study an interplay 
between temperature-induced decoherence and current- 
induced decoherence. For small 8 and V the decoherence 
is solely due to temperature, whereas for larger the volt- 
ages it is mainly current noise. For the 8 = 20 and 8 = 40 
curves in Figure El the 71 decoherence rate is at least at 
the smaller voltages temperature-induced. In contrast, 
the important inter-orbital decoherence rate 712 and also 
72 are not even for smaller voltages. 

If we decrease temperature, the main driving source 
for decoherence is the voltage-induced current noise. The 
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FIG. 7: Difference of the occupation of the two orbitals, (ni — 
71-2), vs. bias voltage V at different orbital splitting 8jT% = 
60, 40, 20, 10, 4, 1 (from top to bottom), temperature T /T% = 
2, and symmetric bare couplings gf x = 0.04, pi 2 A = 0.06. 
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FIG. 8: Differential conductance Gi vs. bias voltage V for 
symmetric initial couplings g AA = 0.04, gi£ = 0.06, tem- 
perature T/Tr = 2, and different orbital splittings 8. Here 

2 

and in the following figures, Go — 2^- denotes the conduc- 
tance quantum. Inset: temperature dependence of the differ- 
ential conductance for = 20 {T/T% = 2: dashed line; 
T/T% = 1: solid line ; T/T% = 0.1: dotted line). 



change with temperature is studied in the inset of Figure 
|SJ showing the sharpening of the Kondoesque peaks. At 
the lowest temperature T/T^ = 0.1, we reached the limit 
where the voltage-induced decoherence is dominating ev- 
erywhere, except for the central peak in Figure |3Ji). 

In Figure we show how the conductances change if 
the couplings of the two dots to the leads are different. In 
part c) and d), the stronger coupled dot (dot 1) is lower 



in energy. The emerging single-dot (spin) Kondo effect 
at V ~ is somewhat suppressed, see Figure c) . In 
contrast in part a) and b), it is the weaker coupled dot 
which is lower in energy (dot 2). Then the spin Kondo 
peak for dot 2 is hardly discernible, see Figure [§ld). 

In the experiments of Ref. 0, the inter-dot Coulomb in- 
teraction U12 is much smaller than the intra-dot Coulomb 
interaction Uu. Hence, the initial inter-dot coupling is 
stronger -g 12 (~ i_ - ±) > - ±), as in 

Figures El and El 

And, as in Figures |S| and El b), the single-dot (spin) 
Kondo peak at V = is much weaker (if at all dis- 
cernible) than the pronounced orbital Kondoesque peaks 
at V » ±6. 



F. Results: Conductance for asymmetric coupling 

Let us now discuss the situation of strongly differ- 
ent (initial/unrenormalized) couplings S^ima ^° ^ ne right 
and the left lead: g^fm^ ^ 9m^m 2 - This is typically the 
case in experiments^ and can in extreme cases give rise 
to a very different physical behavior compared to that of 
the symmetric coupling cases. 

In Ref. 0, this asymmetry is even so strong that one 
of the dots is effectively decoupled from the right lead 
(5^(12) = #2(12) = 0), albeit still coupled to the left lead 
via <?2 L ■ Such a situation is shown in Figure ITU1 Since 
dot 2 is decoupled from the left lead, the orbital Kondo 
effect is now only possible at V = +8, not any longer at 

V = —5. This can be understand by means of Figure 
ITTI At V — +5, we can move an electron from the Fermi 
energy of the left lead into dot 2 and, conserving energy, 
the electron from dot 1 to the Fermi energy of the right 
lead, and vice versa. In the RG equations we hence get 
large couplings, resulting in an orbital Kondoesque peak. 
If the coupling to the leads was symmetric we could have 
similar scattering processes moving the electron from dot 
2 to the right lead's Fermi energy and an electron from 
the left lead's Fermi energy to dot 1 at V = —8. How- 
ever, since g^'^) = ® tms secon d scattering process is not 
possible, as was already argued in Ref. 0. Therefore, in- 
stead of two conductance peaks as in Figures |H1 and |51 we 
see only one at V = +8 in Figure HUl which presents re- 
sults for such strongly asymmetric couplings. The nearly 
indiscernible spin Kondo peak and the missing peak at 

V ~ — 8 qualitatively agree with experiments^ As for 
symmetric coupling, the height of the orbital Kondo peak 
at V = ±<5 decreases with increasing 8. 

An interesting phenomena in the case of asymmetric 
couplings occurs if g^ R — as before but ^> g\ L : As 
shown in Figure 1121 we find a negative differential con- 
ductance which is always present after the orbital Kondo 
peak, i.e., for V > 8. Generally such negative conduc- 
tances occur when the hopping parameters of the Ander- 
son model fulfill \tm\ 3> |tii|, |ifl2 -C \th2\- How can we 
understand this unexpected negative conductance? 
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FIG. 9: Differential conduc- 
tance d vs. bias voltage V 
for symmetric initial couplings 
ff F = 0.04, g$ x ' = 0.01, 
9i2 — 0.06, temperature 
T/T% = 2, and different or- 
bital splittings 5. 



FIG. 10: Part (a) and (c) 
are the same as in Figure |5] 
but now for asymmetric initial 



couplings gi L = 0.04, g RR 
0.01, gh L = 0.01, g% R 
and g& 



0, 

0.12. The Kon- 
doesque peak at V = —5 is 
now missing. Instead of the 
conductance G2 which is zero, 
part (b) and (d) show the dif- 
ference in the occupation of 
the two dots. 
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The basic picture is that -although higher in energy- 
dot 2 is increasingly occupied for V > S. Counterintu- 
itively, it becomes even more occupied than dot 1, see 
the lower part of Figure lT^l The reason for this inversion 
of the occupation (in comparison to the equilibrium oc- 
cupation) can be understood from the sketch in Figure 
1111 The scattering process shown in the Figure is ener- 
getically possible: Because V — S > 0, we can transfer 
an electron from the Fermi level of the left lead to the 
empty states above the Fermi level of the right lead and 
still pay the energy 6 to change the local state from dot 1 
to dot 2. But the reverse process, starting with an elec- 
tron from the right lead's Fermi level, is not possible since 
S — V < 0. All the other processes for changing the local 
state back from dot 2 to dot 1 are much smaller because 
of |izi|/|tfli| <C 1 and |*r2|/|*L2| < 1. Take for example 



the process involving the transfer of the dot 2 electron 
to the right lead and simultaneously that of an electron 
from the left lead to dot 1. The amplitude for this process 
is by a factor of \g^/g^\ (« \t L i/tRi\\t R2 /t L2 \ < 1) 
smaller than that of Figure 10. Hence, the localized state 
is nearly trapped in dot 2, (rii — n 2 ) < 0, for large bias 
voltage. 

Altogether the behavior of the conductance in Figure 
^]can be explained as following: At V ~ S a new, inter- 
orbital channel opens for conductance to the left lead of 
dot 1, the one displayed in Figure ITT1 This explains the 
strong enhancement of G\ for V w 5; the orbital Kondo 
peak in Figure^] With a further increase of V, however, 
there is a dramatic decrease of (ni —712) and hence of the 
number of electrons in dot 1 ((rii)). Because of this, the 
contribution to the current involving solely dot 1, i.e., 
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FIG. 11: Scattering process which connects the Fermi energy 
of the left lead (at +V/2) and dot 2 (upper level) and that 
of the right lead (at — V/2) with dot 1 (lower level); there is 
a separate Fermi level for each dot as indicated. The process 
displayed results, together with the reverse scattering process, 
in a Kondoesque resonance at V ft! 8. For V > S the scattering 
process shown is still possible but not the reverse one. 



the first line in Eq. @, rapidly decreases. This decrease 
dominates for V > 8, resulting in a negative differential 
conductance. 

Negative conductances for quantum dots with many 
levels have been observed experimentally^! and described 
theoretically^ albeit in the sequential tunneling regime. 
Our results show that similar effects are also possible in 
the Kondo regime. 



G. Results: Linear conductance 

Substantial simplifications arise in the linear response 
regime for small voltages since there are no nonequilib- 
rium effects any more like the current induced decoher- 
ence rate 7, and the frequency dependence of the ver- 
tex. This linear conductance can be used as an indica- 
tor whether we have a spin or an orbital Kondo effect, 
corresponding to the limits gi S> <?i2 and 512 3> gi, re- 
spectively. In both cases, the linear conductance for dot 1 
increases at 6 «0, since then dot 1 becomes occupied and 
conducting. But while at larger values of 5 (S/Tg 3> 1) 
the spin Kondo effect still works, the orbital Kondo effect 
breaks down in the linear conductance regime. In Fig- 
ure 1131 a) , the conductance has hence a plateau form in 
the former case, while it has a peak form for the orbital 
Kondo effect. The reason for this dramatic decrease of 
the conductance is simply that the RG flow of the inter- 
orbital coupling (and that for the high energy dot) is 
cut-off at an energy scale ~ 6, preventing a strong renor- 
malization and, hence, suppressing the conductance at 
larger values of 5. 
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FIG. 12: Same as on the right hand side of Figure 1101 but 
now for gi L = 0.01 < g RR = 0.04, and g$ = 0.06. (g£ L = 
0.01, g RR = as before); shown are two different values of 
S and temperature. For these parameters, the differential 
conductance is negative after the Kondoesque peak at V ~ 5 
which can be explained by the inversion of the occupation 
numbers, shown in part (b). 



This different behavior allows us to identify the change 
of the conductance with S, i.e., dG(S)/d6, as a good can- 
didate to determine whether we have a spin or orbital 
Kondo effect. As Figure ED b) shows, this quantity is 
symmetric for the spin Kondo effect but asymmetric for 
the orbital Kondo effect. Note, that dG(8)/dS is often 
accessible experimentally; the two dots need only be con- 
trollable by different gate voltages for changing 6. 



IV. CONCLUSION 

In conclusion, we reported on the extension of the 
nonequilibrium perturbative RG scheme to multi-orbital 
applications. The method keeps the simplicity of the 
poor man's scaling approach, but can be used for 
nonequilibrium transport by including the frequency de- 
pendence of the vertex. However, it is restricted to the 
weak coupling regime. Either current or temperature has 
to cut off the flow to the strong coupling fix point. 

Specifically, we calculated conductances for the case 
of two quantum dots coupled only electrostatically. For 
SU(2) spin symmetry one spin (at V » 0) and two orbital 
Kondoesque peaks (at V « ±5) are found, in agreement 
with the expectations. For the typical experimental sit- 
uation that the intra-dot Coulomb interaction is much 
stronger than the inter-dot Coulomb interaction (trans- 
lating to initial couplings 1712 3> gi, 172) however, the spin 
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FIG. 13: Part (a): Linear conductance Gi(8) as a function of 
orbital splitting 8 for symmetric coupling (see legend; T/T^ — 
2, note that T% depends on g). Part (b) shows dGi{8) / 88 
vs. 8 which is symmetric/antisymmetric for a predominantly 
spin/orbital Kondo effect. 



Kondo peak is much less pronounced, if at all discernible. 

For strongly asymmetric couplings with the extreme 
situation that one of the dots is decoupled from one of the 
two leads, the conductance peak at V ~ — S disappears, 
as observed experimentally. 4 An unexpected result was 
the negative differential conductance immediately after 
the remaining Kondoesque peak at V « 6, occurring for 
initial couplings g RR > g\ L and g RR < g% L . This is a 
genuine nonequilibrium phenomenon. 

Already in the the linear conductance regime, we can 
distinguish between spin and orbital Kondo effect via the 
derivative dGi{5) / 85 which is symmetric and antisym- 
metric w.r.t. 8 — > — <5, respectively. 
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Near resonance u>o, Xm 1 m 2 { u >) behaves like: 

1 

Xm im2 (w) 



(lo - u ) + ry mif 



(A.2) 



where 7 mim2 is just the decoherence rate wanted. 

For the exact diagrammatic evaluation at finite bias 
voltage, one should consider the whole Keldysh contour 
like in Ref. However, for the lowest order contribu- 
tion, it is sufficient to use only the equilibrium Green 
functions as in Ref. [^J Moreover, the proper decoher- 
ence rates to the lowest order in J can be obtained from 
the equation of motion method like in Ref. with a 
truncation of the Eqs. in second order in J. 

The equation of motion of the correlation functions for 
Hamiltonian is written as: 

= (imimi _ r!m2) + (em _ £rai)((imim , im!mi | u 
, V""* tAA' //vroim' + vm 2 mi\\ 
Xk;X>rn'k> 

~ y ] Jmm\ ((-^ 
Xmk\X' k' 



C X'm 1 k l C Xmk', - 



/ tAA' _ tAA' w/ymii 
" / j \ u m-?.m,' J m 1 rn')\\ yv 



'Cw. 



c 

X'm'k' Xm'k' 



X' 



Xk-X'm'k' 



Here, we used the fact that X 



mim2 j^msmA 



VI 



in the single electron subspace. The terms 



on the right hand side can be expressed in a similar way 
by another equation of motion. We then decouple this 
equation of motion by using the approximation 



<(A" lim2 c+C2C+c 4 ;A m2mi }) w p 
[(c24)(ctc 4 ) + (cic 2 )(4c 4 )}((X 



m l m 2 . j£ TTl $171 1 



.(A.3) 



This way, we neglect terms to order J 3 , i.e., our decoher- 
ence rate is only valid to order J 2 . The (cf 02)^04} term 
can be absorbed in a renormalization of the energies 

e m -» e m « e m +y\J^ n f x (e k }i-y j J^fx{e k )\- (A.4) 



Afe 



We can effectively include this renormalization by denot- 
ing with e m not the original level energy of the Hamil- 
tonian but the renormalized e m which is also measured 
in the experiment. The RG Eqs. would include the same 
kind of change if higher order terms were included. 



APPENDIX: DECOHERENCE RATES 

In this Appendix, we calculate the decoherence rates 
via the susceptibility or the corresponding correlation 
function: 

x™H = {{X m ^-X m ^))^ (A.l) 
where ((A; B)) = -i@(t)([A(t), B]) and ((...))„ denotes 
the Fourier transformation to frequency space. 



(u>+e mi -e m2 ) Xmi m 2 (u) = (X m ^-X m ^) 

+inS mim2 J2 C™M«A mm ; A" lim % (A.5) 

m 

plus an additional real part which typically results in a 
modification of e but is of no interest for the decoherence 
rates since these correspond to the imaginary part given 
by 
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B mi m 2 { UJ ) —PO / y {JmimJmmif^( e k)[{ fx'i^k w + fra 2 ^m)] + [J mim {J mim ^ m2m ) + W ^mimi 
W km 

^^-■^J-^^ - w + e m2 -e mi )]} + (mi «-»m 2j w->-w), (A.6) 

(w) =A)^ 1 J lAfe)( 1 - + ^ + f - - + ( w - -w)- 



For mi 7^ m2, we directly obtain the decohercnce rate 
from Eq. i|A.5|) . taking u) — » o; ~ e m2 — ^mi : 

(wo). (A.7) 

This is the rate entering the RG Eq. @ as a cut-off. 
In our perturbative RG calculation, we replace the bare 
couplings pqJ by the renormalized g(co)'s in Eq. (|A.6(I . 
Then, scaling corrections arc included in a similar way 
as in the RG Eqs. themselves. 

For mi = m 2 , the last term in Eq. (|A.5(1 cannot be ne- 
glected. The calculation becomes slightly more difficult 
since a corresponding equation for ((X mm ; X mimi J) ul (the 
terms entering the longitudinal susceptibility) is needed. 
However, often, one does not need the mi = m 2 deco- 
herence rate explicitly. For example, the calculations of 
the paper are for a SU(2) spin symmetry for both lev- 
els/dots. Then, the longitudinal susceptibility and deco- 
herence rate 7 TOlWl with mi = {i,cr} is simply equal to 
the transversal rate 7 miTO2 for m 2 = {i, -a}. 

As an example, the relaxation rate of the SU(N) model 

is 

7 = 2ttN P P 2 T, (A.8) 

which can be obtained by simply substituting J mi m 2 — J 
and J mim2 = J/N. 



An additional check is the spin-i Kondo model. In 
a magnetic field and at a finite bias, we obtain from 
Eq. (|A.5() and an analogous equation of motion for 
(^X mm - X TOimi )) w the transversal and longitudinal spin 
relaxation rates: 

IX = \PV 2 E Jde [h(e - ^)(1 - fy(e + f)) 

XX's z 1 

+fx(e)(l-fy(e))], (A.9) 

7|| = f PP 2 E Jde [/ A (e-^)(l-/v(e+^))]. (A.10) 

This is just the lowest order results one expects from 
nonequilibrium perturbation theory^, i.e.: 

7± n 7|| n Anpl J 2 max{T, ^} (A.ll) 

for max{y,T} > S, and 

7± « |- « Trp^j^j (A.12) 

for max{F,r} < <5. 



1 See, e.g., L. I. Glazman and M. E. Raikh, Pis'ma Zh. Eksp. 
Teor. Fiz. 47, 378 (1988) [Sov. Phys. JETP Lett. 47, 452 
(1988)]; T. K. Ng and P. A. Lee, Phys. Rev. Lett. 61, 1768 
(1988); S. Hershfield, J. H. Davies, and J. W. Wilkins, 
Phys. Rev. Lett. 67, 3720 (1991); Phys. Rev. B 46, 7046 
(1992); J. Konig, J. Schmid, H. Schoeller, and G. Schon, 
Phys. Rev. B 54, 16820 (1996). 

2 Y. Meir, N. S. Wingreen, and P. A. Lee, Phys. Rev. Lett. 
70, 2601 (1993). 

3 D. Goldhaber-Gordon, H. Shtrikman, D. Mahalu, D. 
Abusch-Magder, U. Meirav, and M. A. Kastner, Nature 
391, 156 (1998); S. M. Cronenwett, T. H. Oosterkamp, and 
L. P. Kouwenhoven, Science 281, 540 (1998); J. Schmid, 
J. Weis, K. Eberl, and K. v. Klitzing, Physica B 256-258, 
182 (1998). 

4 U. Wilhelm and J. Weis, Physica E 6, 668 (2000); U. Wil- 
helm, J. Schmid, J. Weis, and K. v. Klitzing, Physica E 9, 
625 (2001); 14, 385 (2002). 

5 A. W. Holleitner, R. H. Blick, A. K. Huettel, K. Eberl, and 



J. P. Kotthaus, Science 297, 70 (2002); A. W. Holleitner, 
A. Chudnovskiy, D. Pfannkuche, K. Eberl, and R. H. Blick 
Phys. Rev. B 70, 075204 (2004). 

6 Also note the case of carbon nanotube quantum dots where 
the band structure implies a two-fold degeneracy, lead- 
ing to multi-orbital Kondo physics as demonstrated by P. 
Jarillo-Herrero, J. Kong, H. S. J. van der Zant, C. Dekker, 
L. P. Kouwenhoven, and S. De Franceschi, Nature 434, 
484 (2005); also see M.-S. Choi, R. Lopez, R. Aguado, 
cond-mat/0411665 

7 J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 
(1966). 

8 B. Coqblin and J. R. Schrieffer, Phys. Rev. 185, 847 
(1969). 

9 In second order scaling theory [one-loop perturbative 
renormalization group (RG)] the prefactor is given by the 
conduction electron bandwidth D of the Coqblin-Schrieffer 
model. It gets corrected by a factor (N poJ)^ 1 ^^ in third 
order scaling theory or when calculating the Anderson im- 



13 



purity model exactly by Bethe ansatz^- 

A. C. Hewson, The Kondo Problem to Heavy Fermions, 
Cambridge University Press (Cambridge, 1993). 
Proposals to artificially arrive at an SU(4) symmetric con- 
figurations include: a triangular quantum dot, see G. 
Zarand A. Brataas, and D. Goldhaber-Gordon, Solid State 
Commun. 126, 463 (2003); the coupling of a small one-level 
dot to a larger grain, see K. Le Hur and P. Simon, Phys. 
Rev. B 67, 201308(R) (2003); and the tuning of the Fock- 
Darwin levels in a vertial quantum dot, see S. Sasaki, S. 
Amaha, N. Asakawa, M. Eto, and S. Tarucha, Phys. Rev. 
Lett. 93, 017 205 (2004). 
12 See, e. g., T. Pohjola, J. Konig, M. M. Salomaa, J. 
Schmid, H. Schoeller, and G. Schon, Europhys. Lett. 40, 
189 (1997); T. Pohjola, H. Schoeller, and G. Schon, Euro- 
phys. Lett. 54, 241 (2001); D. Boese, W. Hofstetter, and 
H. Schoeller, Phys. Rev. B 64, 125309 (2001); W. Izumida, 
O. Sakai, S. Tarucha, Phys. Rev. Lett. 87, 216803 (2001). 

3 W. Izumida, O. Sakai, and Y. Shimizu, J. Phys. Soc. Jap. 
67, 2444 (1998). 

14 D. Boese, W. Hofstetter, H. Schoeller, Phys. Rev. B 66, 
125315 (2002). 

15 L. Borda, G. Zarand, W. Hofstetter, B. I. Halperin, and J. 
von Delft, Phys. R ev. Lett. 90, 026602 (2003). 

16 M. Eto, [cond- mat/0408159] (unpublished). 

17 J. Ye, Phys. Rev. B 56, R489 (1997). 

18 R. Lopez, D. Sanchez, M. Lee, M.-S. Choi, P. Simon, and 
K. Le Hur, Phys. Rev. B 71, 115312 (2005). 

19 A. Rosch, J. Kroha, and P. Wolfle, Phys. Rev. Lett. 87, 
156802 (2001). 

20 J. Paaske, A. Rosch, J. Kroha and P. Wolfle, Phys. Rev. 
B 70 155301 (2004). 

21 A. Rosch, J. Paaske, J. Kroha, and P. Wolfle, Phys. Rev. 
Lett. 90, 076804 (2003); J. Phys. Soc. Jap. 74, 118 (2005); 
A. Rosch, T. A. Costi, J. Paaske, and P. Wolfle, Phys. Rev. 
B 68, 014430 (2003). 

22 Let us also call the reader's attention to another nonequi- 
librium scaling approach, using the flow equation method: 
S. Kehrein, |cond-ma t / 0410341 (unpublished). 



23 P. W. Anderson, J. Phys. C 3, 2436 (1970). 

24 J. Hubbard, Proc. Roy. Soc. A 285, 542 (1965). 

25 L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964) [Sov. 
Phys. JETP 20, 1018 (1965)], J. Rammer and H. Smith, 
Rev. Mod. Phys. 58, 323 (1986). 

26 J. Paaske, A. Rosch, and P. Wolfle, Phys. Rev. B 69, 
155330 (2004). 

27 A difference to Ref. |^ is the way how we deal with the 
frequency conservation at the two vertices of Figure [5] In 
Ref. I2H the conduction electron frequency of the interme- 
diate state is set to the band edge lu' — ±D so that the 
pseudo Fermion line (dashed line in Figure |2J is off reso- 
nance, i.e., its frequency 7^ e m where m is the intermediate 
pseudo Fermion (localized) state. In contrast in Eq. @, 
the localized state is on resonance at e m so that the con- 
duction electron frequency lu' 7^ D. This is similar to the 
traditional poor man's scaling approach where lu' 7^ D and 
the conduction electron momentum is such that et = ±D, 
see, e.g. Ref. [lOJ for details. 

28 This is the Kondo temperature which the two dot system 
had if the orbital splitting was zero but the g's were the 
same. The reason for taking Tk(S = 0) as the energy unit 
is to study changes w.r.t. 5, keeping the other parameters 
constant. 

29 See, e.g., M. Pustilnik and L. I. Glazman, Phys. Rev. B 
64, 045328 (2001). 

30 More details on this issue will be discussed in A. Rosch, J. 
Paaske, J. Kroha, and P. Wolfle (unpublished). 

31 A. T. Johnson, L. P. Kouwenhoven, W. de Jong, N. C. van 
der Vaart, C. J. P. M. Harmans, and C. T. Foxon, Phys. 
Rev. Lett. 69, 1592 (1992); J. Weis, R. J. Haug, K. v. 
Klitzing, and K. Ploog, Phys. Rev. Lett. 71, 4019 (1993). 

32 D. Weinmann, W. Hausler, and B. Kramer, Ann. Physik 
(Leibzig) 5, 652 (1996); M. H. Hettler, H. Schoeller, and 
W. Wenzel, Europhys. Lett. 57, 571 (2002); G. Michalek 
and B. R. Bulka, Eur. Phys. J. B 28, 121 (2002). 

33 W. Gotze and P. Wolfle, J. Low Temp. Phys. 5, 575 (1971). 



